Numerical simulations of fundamental processes in cavity QED: Atomic decay 
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We present results of numerical investigation of a microscopic dynamics of a two-level atom 
embedded into a "linear crystal" of other two-level atoms. These additional atoms play a role of a 
material media. All atoms interact with a multimode cavity field. We study how the decay of the 
initially excited atom is affected by the presence of material media and spectral properties of the 
cavity field. 



I. INTRODUCTION 

Quantum electrodynamics (QED) lies at the heart of 
modern quantum theory. QED is a well established and 
experimentally confirmed theory jl]^] but even fifty years 
after its foundation many features of the atom-field in- 
teraction remain extremely intriguing. In particular, the 
character of the atom-field interaction can be substan- 
tially modified in confined spaces (e.g. the high-Q cavity 
of a micromaser) due to the fact that local properties of 
the electromagnetic (EM) field depend on space bound- 
aries. The radiative properties of atoms and the EM field 
in confined spaces have been investigated for various cav- 
ity QED systems 

Quantum electrodynamics is a local theory, which 
means that the dynamics of atoms and electrons depend 
on local properties of the electromagnetic field with which 
they interact. But local properties of the electromagnetic 
field depend also on conditions imposed by the bound- 
aries of the space region in which the field is confined 
|u| . These conditions are reflected in the quantization 
of the field. We can either quantize the electromagnetic 
field in free space or in a "quantization box" of linear di- 
mension L. Quantum electrodynamics in a box describes 
effects associated with processes inside high-quality (per- 
fect) cavities. In addition, quantization in a box can be 
considered as an approximation to free-space quantiza- 
tion and the two theories must give the same results in 
the limit L — > oo. 

In the first quantum-mechanical description of sponta- 
neous decay of a two-level atom in free space, Weisskopf 
and Wigner [|1J started their calculations with the cav- 
ity modes quantized in a box, and then at a certain stage 
of the calculation, a limit to a continuum of modes was 
taken. This approach gives in first approximation cor- 
rect results (exponential decay of the excited level of the 
atom). The interaction of a two level atom with dis- 
crete (cavity) modes has been described systematically 
by Hamilton |L2| who solved the emission and scatter- 
ing problems exactly for a cubic box by diagonalization 
of the total atom-field Hamiltonian. Later this approach 
was utilized by other researchers for a detailed investiga- 
tion of spontaneous emission of two-level atoms (see for 



instance papers by Davidson and Kozak |13|| and Swain 
|]l4||). Essentially, in all these papers on the spontaneous 
decay of a two-level atom in a cavity (box), the cou- 
pling constant between the atom and the cavity modes 
has been taken to be position independent. This argu- 
ment is perfectly justifiable in free space, when trans- 
lational invariance is valid. On the other hand, when 
the atom interacts with discrete cavity modes in a con- 
fined space the position dependence of the coupling can 
play a significant role. The investigation of this problem 
is not only of theoretical interest. Recent advances in 
experimental techniques have allowed one to study fun- 
damental processes in cavity quantum electrodynamics 
(cavity-QED) |-@,|,| and to verify various effects of 
the atom-field interaction in confined spaces as predicted 
by Schelkunoff Purcell @, Barton and others. 
One of the fundamental processes of cavity QED rep- 
resents the spontaneous decay of a two-level atom. It is 
well known that spontaneous emission from an atom posi- 
tioned very close to the cavity mirror can be significantly 
suppressed. This effect is called the inhibition of sponta- 
neous emission [Q. The deviation from the exponential 
Weisskopf-Wigner decay of the atom in free space |lq | 
has been demonstrated in a number of experiments W. 
Many other interesting questions arise for these QED sys- 
tems. For example, one could ask about the influence of 
cavity mirrors on the dynamics of the atom and the role 
of the position of the atom on the appearance of Poincare 
recurrences (i.e., re-excitations of the atom by radiation 
reflected by the cavity mirrors) (T^| . While the exponen- 
tial character of the decay is not affected by variations of 
the position of the atom around the center of the cavity 
(providing the cavity is large enough), the Poincare re- 
currences depend very sensitively on the position of the 
atom inside the cavity. Namely, variation of the position 
of the atom within a wavelength of the resonant atomic 
frequency can result in an almost complete suppression 
of the first Poincare recurrence of the excited level of the 
atom [^9|. This means that the atom effectively does 
not "feel" wave packets reflected from cavity mirrors for 
times much longer than the time necessary for emitted 
light to "travel" to the mirrors and back to the atom. 
Another example concerns the atom which is positioned 
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close to one of the mirrors. In general one may expect to 
see inhibition of the radiation. Nevertheless, taking into 
account the position dependence of the field-atom inter- 
action it turns out that for some specific distances from 
the mirror (e.g. one quarter of the resonant wavelength 
of the radiation field), the atom decays even faster than 
in free space Q. 

A valuable first insight into modification of the spon- 
taneous emission of the atom into vacuum field is offered 
by the Fermi golden rule M 



employed here. In Section VII we present the convolu 



(1) 



The spontaneous emission rate T a is directly proportional 
to the density of field modes p{oJ a ) in frequency domain 
at the atomic transition frequency uj a \ Vfi is the matrix 
element of the corresponding transition. The presence of 
boundaries (e.g., in the case when the atom is inserted 
into a high-Q microwave cavity) changes the local density 
of field modes and thereby the spontaneous emission can 
be suppressed or enhanced. However, it is by no means 
necessary to change the boundary conditions of the EM 
field in order to modify the spontaneous emission rate. 
This goal can be achieved easily by the presence of other 
atoms which can take part in absorption and re-emission 
of the radiation field. One important example is that of 
an atom embedded in a dielectric host. 

The main task of our investigation here concerns atom- 
field interactions in confined geometries. Starting from 
"first principles" we simulate the dynamics of a system 
of atoms in a cavity. In particular, we consider a cavity 
filled with a "crystal" composed of two-level atoms. We 
investigate the modification of the spontaneous emission 
and the propagation (scattering) of photon wave packets 
within this "crystal". We explore the influence of the 
position dependence of the atom-field coupling on the 
dynamics of the system. In this way also the emission 
and absorption in photonic band gap structures (PBS) 
po| with few atoms can be analyzed. Our micro- 
scopic model based on "first principles" provides us with 
a deeper understanding of the atom-field interaction and 
offers a framework to study systematically the transition 
from microscopic to macroscopic (phenomenological) de- 
scriptions of the systems under consideration. 

In this paper we focus our attention on the modifica- 
tions of spontaneous emission. We describe the model in 
Section O. In Section III we study position dependence 
of the decay of a single two-level atom and in Section |^ 
we discuss the effect of inhibition of spontaneous emis- 
sion. We also analyze how the dynamics of an initially 
excited atom is modified in the presence of other initially 
deexcited atoms in the cavity which play the role of a 
dielectric (see Section [v]) and we study in detail the time 
evolution of atomic populations and quantum-statistical 
properties of the multimode cavity field. In particular, 
in Section VI we focus our attention on the spectrum of 
the field. In addition we will also discuss specific tech- 
nical questions such as the role of the frequency cut-off 



tionless master equation describing the dynamics of the 
initially excited atom in dielectrics interacting with mul- 
timo de c avity field. We summarize our results in Sec- 
tion 
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II. THE MODEL 

We consider a simple one-dimensional model of a cavity 
in which two-level atoms interact with the cavity modes 
in the dipole and the rotating- wave approximations. To 
simplify the model, we neglect all mechanical effects of 
the cavity field on the atom (i.e., the mass of the atom 
is assumed to be infinite). Here we note that this 1-D 
model not only reflects the main features of atom-field 
interaction but also can be mapped onto an isotropic 3-D 
model. 

Under the assumption of perfectly reflecting mirrors, 
the operator of the electric field inside the cavity in the 
Coulomb gauge can be expressed as p|,p3|.p4| 
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E{r) = ^£„e A (a n>A + a^ A J sin(fc n r) 
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where k n — Lo n /c = rnr/L and £ n = ylj^f ■ The two 
orthogonal polarization vectors e\ (A = 1,2) lie in the 
plane perpendicular to the cavity axis; a ny \ and a n < are 
annihilation and creation operators of the n-th mode. 

The Hamiltonian describing the free cavity field can be 
expressed as 
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A n=l 



where we have omitted the zero-point contribution 
^En w ' 1 /^ Summation over discrete modes in Eq.(j|) 
is performed only up to n — N, which means that in our 
model we assume a cutoff for the cavity modes. 

The Hamiltonian describing a set of M non-interacting 
("free") two-level atoms with transition frequencies w^' 
can be expressed as 
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where a 
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|<?) ,-(<?|; \e)j and \g)j denote the 



upper and lower atomic states, respectively. 

When the radius of the atom is much smaller than 
the wavelength of the resonant electromagnetic radiation 
then the atom-field interaction can be described within 
the electric-dipole approximation, i.e., Hint = —d ■ E. 
Neglecting for simplicity all polarization effects, the re- 
sulting interaction Hamiltonian in the rotating-wave ap- 
proximation (FWA) reads 
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where the Pauli spin-flip operators are — \e)j(g\ and 
a_ — \g)j(e\. The position dependent coupling con- 
stants <j4 are given by the expression 

\ 1/2 
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where aig denotes the dipole matrix elements of the 
atoms. The position dependence of the atom-field cou- 
pling constant (^J) given by space-mode functions f n (r) — 
sin(fc n r) may significantly affect the atomic dynamics. 
The total Hamiltonian of the form 
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describes the system of M two-level atoms interacting 
with TV discrete field modes in 1-D cavity. This model 
can be solved exactly because the total number of exci- 
tations 
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is an integral of motion, i.e. [R, H tot ] = 0. 

However, it is impossible in general to find a closed 
analytical solution for the system under consideration. 
Just a few particular cases can be solved analytically, 
and among these is the well known Jaynes-Cummings 
model ^| which describes the dynamics of a two-level 
atom interacting with a single mode cavity field. There- 
fore from the very beginning our treatment will be based 
on a numerical simulation of the cavity QED system. 
Our numerical treatment allows us to calculate proper- 
ties of individual atoms and modes of the EM field, i.e., 
it retains a complete microscopic picture of the problem. 
The Schrodinger equation for the Hamiltonian (^) can be 
transformed into a set of coupled linear differential equa- 
tions for the amplitudes of component states (e.g., basis 
of eigenstates of bare systems) in a finite-dimensional 
subspace of the Hilbert space. Numerical solutions can 
be found using standard methods, e.g. by Runge-Kutta 
methods. Alternatively one can apply the direct diag- 
onalization of the Hamiltonian (0). Numerical solution 
allows us to investigate processes with low initial exci- 
tation numbers (e.g. initially just one or two atoms are 
excited). The number of cavity modes can in this case 
be thousands and the total number of atoms can be up 
to hundreds. 

Using numerical methods, we can analyze the time evo- 
lution of the mean values of the following observablcs: 
i) The occupation of the upper level of the j-th atom 
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ii) The amplitude of the electric field 

N ft, x 1 / 2 
E(r) = J2[^Z) [ a » + 4]sin(fc„r). (10) 

iii) The number of excitations of the fc-th mode 

S(k)=ala k , (11) 

which is used to study the spectrum of the radiation. 

iv) To analyze the space-time propagation of radiation 
wave packets we evaluate mean values of the normally- 
ordered operator for the energy density 



i(r)=:e E 2 (r) 



(12) 



Here normal ordering (the colons above) is adopted to 
eliminate the vacuum-state contribution to the energy 
density of the emitted radiation. 

In what follows we demonstrate the main features of 
the atom-field interaction in confined geometries and de- 
scribe these effects: 

a) Modification of spontaneous emission of the atom in 
the cavity due to the position dependence of the atom- 
field interaction. A partial re-excitation of the atom 
caused by the back reflected radiation (Poincare recur- 
rences) . 

b) Decay in a "crystal" : modification of spontaneous 
emission due to the presence of other atoms, which are 
initially deexcited (the decaying atom can be consid- 
ered as being embedded in a dielectric "crystal" which 
is formed by other atoms). 

c) A model of quantum measurement: a system of two- 
level atoms serves as a device to measure the field spec- 
trum of the considered configuration. The transient char- 
acter of the spectrum can be observed. 



III. DECAY AND RE-EXCITATION OF ATOM 

Within the Weisskopf-Wigner theory Jll| in free space, 
an initially excited atom which is coupled to a contin- 
uum of vacuum field modes decays exponentially to its 
ground state. Representing the usual ID continuum with 
the discrete model (||) for a large cavity (L — > oo), the 
population of the excited atomic level P e decays expo- 
nentially with the rate T a given by the Fermi golden rule 
(HI), i.e., 



P e (t) = exp(-r«t), 



r a = 



oj a \d. 
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In ID "free" space (L — > oo) the decay process is accom- 
panied by the emission of two wave packets (represent- 
ing the one-photon state) propagating to the left and to 
the right from the atom. In the case of the "left-right" 
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symmetry of atomic-wave functions in 1-D (this corre- 
sponds to spherical symmetry in 3-D) each of the two 
wave packets carries half of the initial excitation. This 
process is irreversible as the energy cannot be reabsorbed 
by the atom. In confined geometries the situation differs. 
First, the density of discrete modes is changed due to the 
boundary conditions. The translational symmetry is lost 
and the coupling is position-dependent. In particular, 
when the atom is positioned at the center of the cavity it 
is coupled only to odd modes of the field (for even modes 



the coupling constant ( 
see below). Second, t 



is equal to zero; for more details 
;he two wave packets are reflected 
back by the cavity mirrors and can be (partially) reab- 
sorbed by the atom. This partial restoration of the initial 
state of the atom, the so called Poincare recurrence, can 
be viewed as a consequence of constructive quantum in- 
terference (see below). 

Figure |l| presents the time evolution of the probability 
of the atomic excitation for four different values of the po- 
sition of the atom around the center of the cavity, namely 
An = n - § = 0,±^| ± ^l ; ±^l. To the case when the 
initially excited atom is positioned in the cavity center 
(dotted line) we will further refer as the "free-space" de- 
cay. The central atom interacts only with the odd modes 
and thus the density of (interacting) modes equals to . 
The corresponding "free-space" decay rate r a is given by 
Eq.©. 
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FIG. 1. The time evolution of the population P e (t) of the 
excited atomic level. The atom is shifted from the cavity 
center by Ari = (dotted line), Ari = ±A a /4 (dashed line), 
Ari = ±A a /8 (dashed-dotted line), and Ari = ±A a /16 (solid 
line). The atom is initially prepared in its excited state and 
the multimode cavity field is in the vacuum. The choice of the 
cavity length (in dimensionless units) L = 2ir, the squared 
coupling constant with mode-space function g\ = 1/2 (for 
space-mode function equal to unity) and the atomic transi- 
tion frequency ui a — 100 leads to r a = ir and X a = L/50. The 
first Poincare recurrences appear at the time tE, = 2tt. The 
upper cutoff on frequencies is set to uj cut = 200. 
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From Fig. [T] we see that the first "exponential" stage 
of the decay is (almost) position independent. Providing 
the atom is "far" from the mirrors, i.e., min(ri, L — r\) 3> 
c/r a , the reflected wave packets do not influence the ex- 
ponential decay. More precisely, the density of modes is 
doubled when the atom is shifted from the cavity cen- 
ter. Owing to the position dependence of the atom-field 
coupling (|fj|) the even modes start to interact with the 
atom and the interaction with the odd modes decreases. 
Considering the effective squared interaction constant as 
the average for two neighboring (odd and even) cavity 
modes it decreases to the half value of the squared cou- 
pling constant of the interacting modes for the atom in 
the cavity center. It means that even though the atom 
is shifted from the cavity center the Fermi golden rule 
(Q) with the effective squared interaction constant and 
doubled density of modes leads to the same decay rate 
©■ 

For large enough times, the total excitation energy of 
the atom is transferred to the field, which in turn is effec- 
tively in a one photon (one excitation) state represented 
by two EM wave packets propagating towards the mir- 
rors. For finite cavities at time approximately ^ the 
wave packets are reflected by the mirrors and at tg, ~ — 
they approach the atom, which starts to reabsorb the en- 
ergy from the field. We observe the re-excitation of the 
atom (i.e. the Poincare recurrence can be observed). In 
contrast to the "exponential" stage of the atomic decay, 
Poincare recurrences are very sensitive to small position 
shifts of the atom within a wavelength of the resonant 
atomic transition. In Fig. |l| we clearly see that if the 
atom is positioned at the cavity center (r\ = L/2) then 
at time tu ~ L/c a very well pronounced Poincare re- 
currence of the atomic inversion is seen. One can say 
that at the moment when the Poincare recurrence ap- 
pears the atom "sees" the cavity mirrors [^6|. On the 
other hand, with a small shift of the atom from the cav- 
ity center to Ari = the first atomic recurrence is 
almost completely suppressed. For simplicity we con- 
sider the two emitted wave packets (one to the left and 
one to the right) as monochromatic plane waves (at the 
atomic transition frequency and with the group velocity 
c) . The difference of their geometrical paths is equal to 
4?-. This path difference results in destructive interfer- 
ence due to the accumulated phase difference of n. In 
other words, the atom does not "see" the wave packets 
reflected from the cavity mirrors. Obviously, when the 
two wave packets propagate further, then after the sec- 
ond reflection they accumulate a phase difference of 27T so 
the corresponding Poincare recurrence can then be seen 
(i.e. in this case the atom needs an elapsed time which is 
twice as long compared with the situation when r i = ■§■ 
to "see" the cavity mirrors). 

For the case of an atom positioned at Ari = ±% 
the evolution of the atomic inversion is almost indistin- 
guishable from the case where Ari = 0. At the time 
of appearance of the first Poincare recurrence there is a 
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constructive interference of the wave packets. The trivial 
phase shift 2ir results from the difference of the geomet- 
rical paths which is then equal to A a . 

In the case of an atom positioned at Ar± = ±yg- the 
path difference is equal to A a /4 and the first Poincare 
recurrence is intermediate between the extreme cases 
(Ari = 0, ± g ) considered above. Dephasing of the wave 
packets by ^ results in a partially reduced reabsorption. 
More rigorous analysis should take into account the mul- 
timode structure of the wave packets as an additional 
source of dephasing due to the different (Rabi) frequen- 
cies of the modes. We next note that the second Poincare 
recurrence associated with the second reflection from the 
boundaries starts simultaneously in all cases shown in 
Fig. 0. At the time ~ 2t R the wave packets are merging 
in-phase, i.e. their geometrical paths are equal, which 
results in a partial reexcitation of the atom. 
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FIG. 2. A stroboscopic set of plots describing the 
space-time propagation of the mean energy density of the 
cavity field. We assume the same settings as in Fig. [j] with 
the atom located in the cavity center. We see the two wave 
packets propagating towards the cavity mirrors. These wave 
packets have "sharp" fronts. We note that larger the number 
of modes coupled to the atom the sharper the fronts are. The 
length of the tails of the wave packets depends on the life-time 
of the atom. 

Summarizing this part of the description of the dynam- 
ics of the atom, we can say that while the "exponential" 
character of the decay of an excited atom inside a large 
cavity is not influenced by small shifts of the atomic posi- 
tion, the first Poincare recurrence is a position-dependent 
interference effect. The basic features can be explained 
using a simple classical point of view based on two inter- 
fering monochromatic waves. This fact may be thought 



of as rather surprising, as mathematically the Poincare 
recurrences can be related to the discrete nature of the 
cavity modes (with equal frequency spacing) . To be more 
precise, the phase-matching conditions necessary for the 
appearance of Poincare recurrences in the atom-field sys- 
tem can be associated with the evolving phase factors 
e -ttE k f con t r ibuting eigenstates of the total Hamil- 
tonian (Q): a Poincare recurrence, i.e., a partial restora- 
tion of the initial state, can appear at time tn such that 
the relation E^ta ~ 2ir is valid for many values of k (for 
more details see (27]]). 

In Fig. U we present a stroboscopic set of plots describ- 
ing the space-time propagation of the energy density of 
the cavity field for the same physical situation (initial 
state) as considered in Fig. |] when the atom is in the 
center of the cavity. We see two distinct wave packets 
propagating to the right and to the left. Reflection of 
the wave packets from the cavity mirrors (at time t ~ 
is nicely demonstrated and the subsequent re-excitation 
of the atom is synchronized with the interference of the 
wave packets in the center of the cavity (compare with 
Fig. ff). 



IV. INHIBITION OF SPONTANEOUS EMISSION 

In the previous section we have considered situations 
when the atom is "far" from the cavity mirrors [i.e. 
min(ri,L — r\) 3> c/T a ] and the wave packets reflected 
by the mirrors do not directly affect the initial sponta- 
neous decay of the atom. On the other hand, for dis- 
tances between the atom and one of the cavity mirrors 
smaller than c/T a (here l/r a is the spontaneous decay 
time in a free space) deviations from exponential decay 
should be expected In particular, the decay of a 

two-level atom which is positioned very close to the cav- 
ity mirror can be significantly suppressed. This effect is 
called the inhibition of spontaneous emission . The in- 
hibition of spontaneous emission is a position-dependent 
effect which is related to the position dependence of the 
atom- field coupling constant (||). In Fig. || we present 
numerical simulations for the time evolution of the pop- 
ulation of the upper level of the atom described by the 
model interaction Hamiltonian ([5]) . The atom is assumed 
to be initially in its excited state and the field in the 
vacuum. We consider several typical physical configura- 
tions. Firstly, for reference, we plot the atomic decay of 
the atom positioned at r± = X a /8 (solid line) which is 
indistinguishable from the exponential decay of the atom 
at the cavity center (i.e., P e ss e _ra * for t < tn). For 
other atomic positions t\ = A a /16 and r\ — A a /32 (here 
A a = L/50) we clearly see that the closer the atom is 
to the mirror the slower the spontaneous decay is, that 
is the inhibition of spontaneous radiation is transparent 
for the considered positions of the atom. On the other 
hand for very specific atomic positions close to the mir- 
ror the opposite effect - the enhancement of spontaneous 
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emission - takes place. Namely, for the atomic distance 
Ti = A a /4 the atom decays as P e « e~ 2Tat , i.e., it radi- 
ates twice as fast compared with the free-space case [see 
the reference case ri = A a /8]. 
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FIG. 3. The time evolution of the population P e (t) of 
the excited atomic level for the atom very close to the cav- 
ity mirror. The atom is considered at the following posi- 
tions: ri = A a /2 (dotted line), r\ = A a /4 (dashed line), 
ri = A a /8 (solid line), r\ — A a /16 (dashed-dotted line) and 
T\ — Aa/32 (dotted line). The "reference" exponential decay 
of the atom at the cavity center r\ — L/2 coincides with the 
case ri = A a /8. The initial conditions and other parameters 
are as in Fig. |l[ 

The origin of inhibition or enhancement of spontaneous 
emission in the context of the model used in this paper 
relies on the position dependence of the atom-field cou- 
pling (^). In particular, for r\ — A a /4 the spatial-mode 
function sin(fc„ri) ss I for all the modes close to the 
resonant frequency uj a irrespective of whether n is even 
or odd. This means that the density of modes is in- 
creased by a factor of two compared with the case of the 
atom at the cavity center r% — L/2 when the modes with 
even n are decoupled from the atom [sin(fc ra L/2) = for 
even n and sin(fc„L/2) = I for odd n]. The increased 
density of modes implies an enhancement of the spon- 
taneous emission. In a similar way, when r\ = A a /8 
the spatial-mode function sin(fc„ri) ss I/a/2 for all n 
around the atomic transition frequency. However, the 
decrease of the squared interaction constants is compen- 
sated for by an increase in the density of interacting 
modes (compared with the atom at the cavity center) 
and thus the spontaneous emission rate retains the value 
r a . For the other extreme case r\ = X a /2 all modes 
around u a are significantly decoupled from the atom [now 
sin(A:„ri) ss 0] which results in the dramatic inhibition of 
the spontaneous emission as seen in Fig. pi For other 
positions shown in Fig. || (ri = A Q /16,A Q /32) the slow- 
ing of the spontaneous emission rate is given by the fac- 
tor r(ri)/r a m 1 — cos(2fc a ri). We note that in the 



case n = A Q /I6 the atom decays completely while for 
T\ = A Q /32 the exponential decay law is interrupted by 
the Poincare recurrence at 2tR. 
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FIG. 4. The time evolution of the population P e (t) of the 
excited atomic level for the atom at the following positions: 
ri — X a (dashed line); ri = \ a + A a /4 (dotted line); and 
7"i = A a + A a /8 (dashed-dotted line) which are compared with 
the exponential decay of the atom in the cavity center (solid 
line). Other settings are as in Fig. |f|. The suppression and 
the stimulation of the emission caused by the reflected wave 
packet are clearly see. 



For distances of the atom from the mirror larger than 
the wavelength of the resonant transition A a the interfer- 
ence with the reflected wave packet can either stimulate 
or suppress the emission of the atom. To be specific we 
show in Fig. [| the time evolution of the population of 
the upper level of the atom which is positioned at three 
distances n = A a (dashed line); ri = A a + A Q /4 (dot- 
ted line); and r\ = X a + X a /8 (dashed-dotted line) which 
are compared with the exponential decay of the atom in 
the cavity center (solid line) . The phase accumulated by 
the wave packet during the round trip from the atom to 
the neighboring mirror and back is in the case r\ = X a 
equal approximately to 5ir (here the additional contri- 
bution of 7r is due to the reflection from the mirror), 
i.e. there is a destructive interference between the wave 
packet and the atom which results in the suppression of 
the radiation. On the other hand, when r\ = X a + A a /4 
the accumulated phase is approximately 67r, which leads 
to constructive interference. It means that the reflected 
wave packet, when it arrives at the position of the atom, 
starts to stimulate the atomic emission. In the units used 
in this simulation, the arrival time of the reflected wave 
packet is at approximately t ~ 0.16 which coincides with 
the change of the initial exponential decay of the atom. 
When the atom is at the position r\ = X a + A Q /8, the 
accumulated phase of the reflected wave packet is ll7r/2 
which gives rise to a partial suppression of radiation. 
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V. SPONTANEOUS EMISSION IN A "CRYSTAL" 

Atomic radiation can be crucially modified by the pres- 
ence of other atoms in the cavity. Obviously, if the dis- 
tance between the atoms is large enough then the expo- 
nential decay of the originally excited atom is not affected 
significantly. On the other hand when the atoms are 
placed close together the situation is different (one of the 
consequences is a collective behavior of the atoms which 
might result in superradiance, see for instance p8|). 

In this section we consider a specific initial condition 
when the initially excited atom is surrounded by a col- 
lection of two-level atoms in the ground state. These 
additional atoms can be considered as a linear "crys- 
tal" . By changing the density of the atoms we can model 
systems such as atomic structures embedded in optical 
lattices (for interatomic distances comparable with the 
wavelength of the atomic transition) or dielectrics (for 
much smaller interatomic distances). The cavity QED 
system with trapped atoms [29| l represents a quite new 
experimental system. The transfer of excitations between 
particular atoms which are captured in an optical poten- 
tial is mediated by the cavity field [ p0[ . 
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FIG. 5. The modification of the spontaneous emission of 
the atom at the cavity center being surrounded by identical 
atoms which form a linear "crystal" . The central atom is ini- 
tially excited and the others de-excited. The vacuum field 
is in the vacuum. The regular linear crystal which fills the 
central part of the cavity is composed of M = 101 atoms 
with the interatomic distance a — A a /2 (long dashed line), 
a = A a /4 (short dashed line), a — X a /8 (solid line), a = A a /16 
(dot -dashed line). In the case a = (dotted line) all the 
atoms are positioned at the cavity center. Single atom decay 
M = 1 (sparse dotted line) is shown for reference. 

The modification of the spontaneous emission of the 
atom embedded in a linear "crystal" of two-level atoms 
is shown in Fig. ||. The regular crystal built of M = 101 
atoms fills the central part of the cavity. Initially the only 
excited atom is located at the cavity center and decays 



in the vacuum field. The modification of the sponta- 
neous emission depends dramatically on the interatomic 
distance a. In the case of the "lattice" constant being 
a = A a /2 (long dashed line) we can observe strong sup- 
pression of the spontaneous emission while for a = A a /4 
(short dashed line) an enhancement of radiation com- 
pared with the single atom system (dotted line) takes 
place. The origin of this behavior is related to either 
destructive or constructive interference effects, respec- 
tively. From other examples, for a — A a /8 (solid line) and 
a = A a /16 (dot-dashed line) it is seen that by increas- 
ing the density of the linear "crystal" the atom radiates 
more slowly. Moreover, the de-excitation is incomplete, 
as an increasing part of the excitation is captured by 
the initial state. This subradiant behavior has been al- 
ready analyzed for the extreme case when all the atoms 
are located at the same position (e.g., the cavity center) 
pl[ . The initial excitation is captured in the asymmetric 
atomic state and only a small part ~ -j? is radiated in 
the cavity field. 
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FIG. 6. The time evolution of the population P e (t) of the 
excited atomic level for the atom surrounded by identical 
atoms which form a random linear structure. Averages over 
100 random configurations with one atom within the lattice 
constant a = (A a /8) (triangles), a — (A a /4) (circles), and 
a = (A a /2) (squares) are compared with the corresponding 
regular crystals for a = A a /8 (dotted line), a = A a /4 (dashed 
line), and a — A a /2 (dot-dashed line). Setting are as in Fig. ^. 
Single atom decay (M = 1) in "free-space" (solid line) is 
shown for reference. 

The regular "crystal" represents a rather specific and 
idealized case. Positions of atoms can fluctuate due to 
various reasons (for example in the case of optical lat- 
tices with shallow wells formed from optical potentials). 
To simulate the situation when the atoms are not regu- 
larly distributed in the cavity we next consider random 
configurations of the atoms. Specifically, the atoms are 
placed randomly such that within each lattice constant 
there is just one atom. Depending on the particular posi- 
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tions of the atoms, the dynamics of the originally excited 
atom can change dramatically. The atomic radiation can 
be either enhanced or suppressed. To obtain some effec- 
tive "macroscopic" picture from our simulations, we have 
averaged over many random configurations. 

We present the results in Fig. |[ The dashed (dot- 
ted) line in this figure shows the time evolution of the 
atomic population of the initially excited atom when the 
atoms are regularly positioned with the lattice constant 
a = A a /4 (a = A a /8), representing enhancement (sup- 
pression) of radiation with respect to the "free-space" 
decay (solid line). The results of numerical simulations 
corresponding to averaging over many (100) random con- 
figurations of atoms are presented for the average dis- 
tance between atoms (a) = A Q /4 (o) and (a) = A a /8 
(A). In both cases the radiation of the atom is sup- 
pressed compared with the "free-space" decay. Another 
common feature of the dynamics in this case is that in 
both cases the atom does not radiate away completely 
the initial excitation energy. Both effects (suppression of 
radiation and "excitation trapping") are caused by the 
collective influence of the crystal atoms. 

An increase of the density of the linear "crystal" (e.g. 
for a = A a /16) rapidly diminishes differences between the 
regular crystal and the corresponding "random" crystal 
with one atom randomly positioned within the lattice 
constant. For completeness we included in Figjq also 
the case of the regular "crystal" with a = A Q /2 (q) and 
the average over random configurations with (a) = A a /2. 
Here the destructive interference effect which leads to the 
strong inhibition of the radiation for the regular crystal 
is deteriorated for random atomic configurations shown 
in Fig.§ 
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FIG. 7. The total excitation of the atoms Ratoms = ^2 P& 
which form a regular linear "crystal" as in Fig. |H| For the in- 
teratomic distance a = A a /4 also small "crystals" made of 
M = 11 (dotted line) and M = 21 (dot-dashed line) are 
considered. 

We note that the modification of the spontaneous emis- 



sion is a local effect, i.e., the atomic decay is influ- 
enced only by neighboring atoms. To check this we have 
performed simulations with only 10 neighboring atoms 
(M = 11). We have found that the "exponential" stage 
of the spontaneous emission is unchanged comparing with 
the case M — 101, and differences occur only on a much 
longer time scale. 

It is instructive to note that part of the initial excita- 
tion energy is captured in the atomic system. In partic- 
ular, the sum of atomic excitations R a toms = Y] Pe (t) 
depends mainly on the number of atoms M. This ten- 
dency is confirmed in Fig. ^ which shows the cases 
M = 11,21,101 for a = A a /4. In general, the oscilla- 
tion patterns reflect complex interference effects. Never- 
theless one can trace a very general tendency in the pic- 
ture - the "crystal" atoms which surround the initially 
excited atom play the role of semi-transparent mirrors 
placed very close to the atom. Therefore the results par- 
tially resemble the case of the single atom in the vicinity 
of a mirror (compare Fig. [| with Fig. ||). 



VI. SPECTRUM OF THE CAVITY FIELD 

Within the framework of cavity QED when the field in- 
teracting with the atoms is confined within ideal mirrors, 
there is nothing like a stationary regime which is neces- 
sary for the derivation of a time-independent spectrum of 
the field. The spectrum is intrinsically time dependent. 
In this case an operational definition of time-dependent 
spectrum can be given by excitation probabilities of the 
cavity modes [see Eq.(pT[)]. 

The spectrum of the cavity field is affected by the po- 
sition of the atom. In particular, if the atom is located in 
the cavity center even modes are completely decoupled 
from the atom and only odd modes can become excited 
[see Eq.(^)] establishing in this way oscillations in the 
spectrum of modes. However, at the point when the to- 
tal excitation energy of the atom is transferred to the 
field, the envelope of the spectrum is "Lorentzian" irre- 
spective of the position of the atom (providing that the 
decay is exponential). 

On the other hand, it should be stressed that the spec- 
trum of the interacting modes is highly transient even 
during the exponential decay period. It undergoes a 
gradual narrowing from a broad flat spectrum (initially 
all modes are in the vacuum state with the same prob- 
ability) towards a Lorentzian- like line of width r Q . The 
narrowing is accompanied by transient oscillations of the 
spectral envelope. This transient behavior is illustrated 
in Fig. |§| which shows the spectral envelopes at different 
time moments during the exponential decay of the atomic 
excitation. At the time t~2 the envelope of the cavity- 
field spectrum reaches its quasi-stationary shape, being 
very close to the corresponding (Lorentzian) emission 
spectrum usually associated with the free-space emission 
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FIG. 8. Populations of cavity modes (i.e. the spectrum) at times t = 0.3, 0.7, 1, 3 (in units of l/r a ) for the atom located at 
the cavity center. The atom is initially prepared in its excited state and the cavity field is in the vacuum (for other conditions 
see Fig. |l|). For comparison purposes we show the overlaps of the eigenstates of the total Hamiltonian ([?]) with the initial state 
(these overlaps are denoted by the symbol A). The ordering of eigenstates is given by their eigenvalues on frequency axis. The 
even modes which do not interact with the atom have no overlap with eigenstates of the total Hamiltonian (see oscillations in 
the function \{ijj \<f> k )\ 2 ). 



It is worth noticing that there is a close relation be- 
tween the emission spectrum and the "spectrum" of 
squared scalar products (overlaps) between eigenvectors 
|$fc) of the total Hamiltonian (^) and the given initial 
state |^(io)), i-e., 

S e (fc) = |MM|$ fe )| 2 . (14) 

From Fig. || it is evident that the "spectrum" of overlaps 
(shown as •) resembles the emission spectrum of the com- 
pletely deexcited atom. In other words, the "spectrum" 
of overlaps offers an important time-independent charac- 
teristization of the system under consideration. If there 
exists a quasi-stationary spectrum of the cavity modes, it 
should mimic the "spectrum" of overlaps. In addition, a 
shift of the atomic transition frequency in the spectrum 
of eigenvalues can be associated with the energy shift. In 
our calculations we use a broadband approximation, in 
which the frequency dependence of the atom-field cou- 
pling can be neglected. Thus for the symmetrical upper 
frequency cutoff uj cut — 2oj a the energy shift equals zero. 

We next turn our attention to the fact that an addi- 
tional system of two-level atoms with different transition 
frequencies inside the cavity can serve as a device 
to measure the spectrum of the cavity modes. Within 
this model of measurement the initially deexcited (an- 
alyzer) atoms are placed far enough from the decaying 



atom. The coupling of the analyzer atoms to the cav- 
ity modes is so weak that the dynamics of the cavity 
field is essentially unaffected on the relevant time scale. 
It means that the linewidths of the analyzer atoms are 
much smaller than the linewidth r a of the decaying atom. 
The analyzer atoms with different transition frequencies 
(7) 

u>a thus function as narrow frequency filters interacting 
effectively only with the cavity modes on exact resonance 
with particular uia . The upper level excitation proba- 
bilities of the analyzer atoms are proportional to the in- 
tensity of the cavity field at the position of the analyzer 
atoms |32|. In order to map the excitations of analyzer 
atoms to the spectrum of the cavity modes it is necessary 
to demand that the coupling constants (linewidths) of the 
analyzer atoms are equal. In general, the analyzer atoms 
give the local and time dependent frequency spectrum 
which is not necessarily the same as excitations of cavity 
modes ( |ll] ) which give the "instantaneous" spectrum in 
the whole cavity. In the case of the atomic decay the 
two spectra agree if the times they are read out from the 
analyzer atoms and the modes are chosen appropriately. 

The normalized absorption spectrum read out from the 
analyzer atoms is show in Fig. ^|. Here we have considered 
that the decay rates (i.e. linewidths) T of the analyzer 
atoms are mutually equal and T = 10 _4 r a . A bunch of 
one hundred of analyzer atoms is positioned Ar = 0.5 







apart of the decaying atom in the cavity center. The 

(7) 

atomic frequencies uia are equally spaced around the 
transition frequency Lu a of the central atom. We see very 
good agreement with the envelope of the spectrum of the 
cavity modes. More precisely, in Fig. |^ we compare the 
spectrum of the cavity modes at times t = 0.3 and t = 2 
with the absorption spectrum of the analyzer atoms at 
delayed times t + tf where tf is the time of flight of the 
light from the decaying (central) atom to the analyzer 
atoms. 



Spectrum of the decaying atom 




" *ioo* * 
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FIG. 9. The spectrum of the cavity modes at times t = 0.3 
(★) and t = 2 (o) is compared with the normalized absorption 
spectrum of the analyzer atoms at times t — 0.3 + tf (dashed 
line) and t = 2 + tf (solid line). Here tf is time of flight of 
the light from the decaying atom at the cavity center with T a 
as in Fig. [I] to the analyzer atoms with mutually equal decay 
rates F = 10 _4 F a . The distance of the analyzer atoms from 
decaying atoms is Ar = 0.5 and the time of flight tf — 0.5 in 
choosen units (c = 1). 



A. Numerics and cutoff dependence 

Our numerical calculations have been performed in the 
broadband approximation for the interaction constants 
given by Eq.(p), i.e., we have neglected the frequency de- 
pendence of the coupling constants replacing u> n by u> a . 
This approximation is valid only for large enough cavi- 
ties with L » ) a and "weak" interaction regimes with 

A rather subtle point is the choice of the frequency 
cutoff. Strictly speaking, the model interaction Hamil- 
tonian (|5|) with the interaction constants (|^) within the 
broadband approximation leads in second-order pertur- 
bation theory to logarithmically divergent energy shifts 
for w cut — > 00 p4j . (Note that in our numerical cal- 
culations we have eliminated the shift of the excited 
atomic level by choosing a symmetrical frequency cut- 
off oj cu t — 2oj a .) Without the broadband approximation, 



when the frequency dependence of the interaction con- 
stants ([|) is taken into account, the energy shifts diverge 
linearly. It is well known Jl],|| that if instead of the dipole 

approximation Hi nt — —d-E we start with Hi nt — —ft- A 
then after the RWA is applied one obtains the interac- 
tion Hamiltonian (||) but with a different frequency de- 
pendence of the interaction constant, i.e. 




\ heoL J 



(15) 



In the broadband approximation the interaction con- 
stants given by Eq.(|) and Eq. (|l5|) are identical and the 
results do not depend on the choice of the interaction 
Hamiltonian. On the other hand, without the broad- 
band approximation the results are biased by the choice 
of the frequency dependence of the atom-field coupling. 
From the mathematical point of view, the coupling given 
by the expression ( p"5| ) does not lead in second-order per- 
turbation theory to divergent energy shifts for oj cu t — > 00. 
Obviously at the point when the two effective Hamiltoni- 
ans considered above lead to different results one has to 
be careful whether the model is physically relevant (for 
more details see Ref. 0). 



VII. MASTER EQUATION FOR THE ATOM IN 
DIELECTRICS 

The system of atoms and field modes under consid- 
eration in an ideal cavity represents a closed system 
with unitary dynamics governed by Schrodinger equa- 
tion. In this Section we will consider the decaying atom 
as an open system in the environment represented by field 
modes and other initially unexcited atoms. This analysis 
can give us a deeper insight into the problem of dynam- 
ical evolution of the atom in dielectrics interacting with 
many cavity modes. In general, an open system S (in our 
case the atom which is initially excited) interacts with an 
environment E (the other atoms surrounding the origi- 
nally excited atom and cavity modes) Q . In this section 
we consider the archetypal system + environment model 
which is specified as follows: Let 7i s denotes a Hilbert 
space of the system S, and TC E is the Hilbert space 
associated with the environment E. The Hamiltonian 
H SE — H s ®l E +H int + l s ®H E of the composite system 
S(BE acts on H s ®Ti. E . It is assumed that S®E is a closed 
finite- dimensional system which evolves unitarily. The 
density operator p SE (t) of this composite system is gov- 
erned by the von Neumann equation with the formal solu- 
tion p SM (t) = exp[-i(t-to)H SE ]p SE (t )exp[i(t-t )H SE ], 
where the initial state is p SE (to) = P s (io) ® P E (io) an d 
h = 1. The reduced dynamics of the system S is then 
defined as 

p s (t) f(t, t ) Ps (t ) = Tr E [p SE (t)] . (16) 

By definition, T(t, to) is a linear map which transforms 
the input state p s (to) onto the output state p s (t). In our 
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recent paper pm we have addressed the question how to 
determine (reconstruct) the master equation which gov- 
erns the time evolution of the reduced density operator 
p s (t) . It has been shown that this master equation can 
be written in the convolutionless form (we omit the sub- 
script S) 



-p(t) = £(t,t Q )p(t)- 



(17) 



which is possible due to the fact that in the finite- 
dimensional Hilbert spaces matrix elements of density 
operators are analytic functions. Consequently, T(t, to) 
are non-singular operators (except perhaps for a set of 
isolated values of t) in which case the inverse operators 
T(i, to) exist and the Liouvillian superoperator can be 
expressed as 



C(t,to) := 



dt 



T(Mo) 



T- l (t,t ). 



(18) 



We note that T(t,to) is uniquely specified by H SE 
by the initial state p E (to) of the environment. 



and 
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FIG. 10. The time evolution of the time-dependent decay 
rate T(t). We assume the initially excited atom to be in the 
center of the cavity, other settings are as in Fig. Jjj. In the 
chosen units we obtain from the Fermi golden role (see |l3| ) 
the decay rate T a = n (see dashed line). We consider two 
cases: firstly when there is just a single excited atom in the 
cavity (dotted line) and secondly, when the excited atom is 
surrounded by 100 atoms (solid line) which create a "linear" 
crystal with regular spacing between atoms (a — A a /8). 

In an earlier paper [Q we have propose a general al- 
gorithm how to reconstruct the Liouvillian superopera- 
tor C(t, to) from the knowledge of the unitary evolution 
of the composite S © E system [see Eq.(|l6|)]. From this 
knowledge the master equation (17) can be uniquely de- 
termined. The dynamics of the open system (in our case 
the atom) is given exclusively in terms of the system oper- 
ators. Environmental degrees of freedom are completely 



eliminated from the reduced dynamics. Nevertheless, the 
state of the environment may change during the time evo- 
lution due to the interaction with the system. That is, 
there is no need to employ the assumption that the envi- 
ronment is a "big" reservoir which does not change under 
the action of the system. 

To apply the formalism presented in Ref. |34| we re- 
mind ourselves that the initial state of the atom-field 
system considered in the present paper reads 



l*ft>)> = |e)i|fl> 7 |0>. 



fe' 



(19) 



where |e)i describes the excited state of the initially ex- 



cited atom, while \g)-* := |g); 



® \q)m describes the 



rest of the M atoms which are initially in the ground 
state. The vector |0)g denotes the vacuum of the cavity 
field. Because the model Hamiltonian H tot (Q) is chosen 
so that the number of excitations in the system is an in- 
tegral of motion we can express the state vector of the 
atom field system at time t as 



|*(*)> = d(i)|e)i| 5 >j|0>g 

M 

+ £e*(t)l0>i|e J ->j|O>j ? 
+ I>fc(*)l<7>i|<?)j|lfc>S 



(20) 



where \ej)j describes the state vector of a set of M — 1 
atoms out of which the j-th atom is excited, while 
describes the state of the cavity field with the fc-th mode 
in the Fock state |1) and all other modes in the vacuum 
state. Using this notation we can express the master 
equation for the originally excited atom as [B4| 



d , .s(t) r . . , 



(21) 



r(t) 



2cr_/5(7 + — <t + <5-_/5 — pa+ <7_ 



where the time-dependent decay rate T(t) and the time- 
dependent dynamical energy shift S(t) can be expressed 
through the probability amplitude c\ (t) as 



r(i) = Re [r)(t)] ; 5(t) = Im fo(t)] 



where 



v(t) 



1 da(t) 



ci(t) dt 



(22) 



(23) 



In general the parameter rj cannot be derived in an an- 
alytical form. In Fig. [l(] we present results of numerical 
evaluation. We assume the initially excited atom to be 
in the center of the cavity. In the chosen units we obtain 
from the Fermi golden rule [see Eq.(|l^) the decay rate 
r a = 7r (see dashed line). We consider two cases: firstly 
when there is just a single excited atom in the cavity 
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(dotted line) and secondly, when the excited atom is sur- 
rounded by 100 atoms (solid line) which create a "linear" 
crystal with regular spacing between atoms (a = A a /8). 
In the case of the single atom T(t) oscillates around the 
value r a . The amplitudes of these oscillations are rela- 
tively small till the recurrence time when it eventually 
takes negative values (i.e. the atom starts to absorb en- 
ergy from the field - compare with Fig. ^). In the second 
case which mimics the decays of the atom in dielectrics 
the time evolution of T(t) is more complex. At the ini- 
tial instants T(t) oscillates around the value T a , but then 
it rapidly decreases and takes negative values - this is 
correlated with the absorption of energy from the wave 
packets reflected by surrounding atoms (see Fig. ^). 
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VIII. CONCLUSIONS 



In this paper we have numerically studied a micro- 
scopic model of the cavity QED describing an atom in- 
teracting with multimode electromagnetic field. The ini- 
tially excited atom is surrounded by a set of other atoms 
which represent a dielectric "crystal" on a microscopic 
level. We have shown how the decay of the atom is mod- 
ified due to the mode structure of the cavity field and the 
presence of neighboring atoms. 

In the cavity QED model considered here we have ne- 
glected all mechanical effects of the EM field on atoms. 
To account for effects of quantized center-of-mass mo- 
tion one can generalize the model to a system of trapped 
atoms interacting with EM field in the cavity. We can 
assume a situation when due to laser cooling the atoms 
embedded in a cavity (or PBS) are effectively trapped 
by harmonic (optical) potential. There is a close analogy 
with the model of trapped ions (29) . 

In our next paper we focus on the propagation of one- 
photon wave packets in a cavity. These wave packets are 
represented as superpositions of many cavity modes. For 
example, we will study in detail scattering of a photon 
wave packet on a single two-level atom. We will show 
that the atom within this framework effectively behaves 
as quantum beam splitter or semi-transparent mirror. We 
will also discuss dynamics of the wave packet incident on 
a "crystal" composed of two-level atoms which fill some 
part of the resonator and thereby the transmission of the 
energy of EM field depends on the density of the atoms. 
This model will help us to understand on the microscopic 
level how photon wave packets propagate through dielec- 
tric media and to estimate (at least qualitatively) the 
speed of their propagation from first principles. In ad- 
dition the model will allow us to formulate the quantum 
version of the extinction theorem |?5| which on the classi- 
cal level ||(| explains how electromagnetic waves packets 
propagate in material media. 
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